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We know that the length of minor great circle arc joining any two arbitrary points on a sphere of finite radius is 
the minimum distance between those points. Here we are interested in finding out the minimum distance or 
great circle distance between any two arbitrary points on a spherical surface of finite radius (like globe) for the 
given values of latitudes & longitudes. 


Let there any two arbitrary points A(@, ,A1) & B(@z , Az) 
on the surface of sphere of radius R & centre at the point 
O. The angles of latitude @, & m2 are measured from the 
equator plane & the angles of longitude A, &A, are 
measured from a reference plane OPQ in the anticlockwise 
direction. Here, we are to find out the length of great circle 
arc AB joining the given points A & B. Draw the great circle 
arcs passing through the points A & B which intersect each 
other at the peak (pole) point P & intersect the equatorial 
line orthogonally at the points D & C respectively. (As 
shown by the dotted arcs PD & PC in the figure 1) 


Join the points A, B, C & D by the dotted straight lines 
through the interior of sphere to get a plane quadrilateral 
ABCD. 


Now, the angle between the great circle arcs AD & DC, 





subtending the angles a= ¢@, & B =A, —A, respectively 


at the centre O of the sphere, meeting each other at the 
common end point D, is 6 = TI , hence the (plane) angle 
_< ADC = @p between the chords AD & DC of great circle 
arcs AD & DC is given by HCR’s cosine formula as follows 


i ain! ine eee” eas’ cask 
cos0, = sin> sin? + cos>cos > cos 


Figure 1: The dotted great circle arcs passing through two given 
points A & B, intersecting each other at the peak (pole) point P, meet 
the equatorial line orthogonally at the points D & C respectively on a 
spherical surface of finite radius R 


V 0<a,fp,0O<na> 0<6,<84 


Now, setting the corresponding values in the above formula, we get 











cos@y = sin-—sin (“ = “) + cos —cos am “) wos 
. 2 2 2 2 2 
A,—A 1 
cosOy = sin! sin 7 +) V 0<6@)< 2 ieeieieaedes (CL) 


Similarly, setting the values of anglesa=@, & B=A,—-1A, & A= T/> in the formula, the (plane) angle 
_< BCD = 6, between the chords BC & CD of great circle arcs AD & DC, meeting each other orthogonally at 


the common end point C, is given as 


Copyright© H.C. Rajpoot 








_. ba, (“ zy “) p2 (* = *) 1 
cos0, = sin—sin + cos —cos cos > 


2 2 2 2 
A, —A ie 
cos0>; = sin &? sin (254) V0<6@< 2 uieivwi ween (Z) 


Now, join the points A, B, C & D to the centre O of the sphere by the dotted straight lines (As shown in the 
figure 1 above). 


In isosceles AAOD, 


OA=OD=R, —<AOD= 4, 


ot 2 #1 
sin 7 R => AD=2Rsin 2 


Similarly, in isosceles ACOD, 
. (A2— 41 
OC =OD=R, COD=A,-A, > CD= 2Rsin (2) 
Similarly, in isosceles ABOC, 
OB =OC=R, <BOC=¢, > BC= 2Rsin © 


Now, join the vertices B & D in the plane quadrilateral ABCD to get ABCD (As shown in 
the figure 2). Applying cosine rule in ABCD as follows 


BC? + CD? — BD? 


ABCD = 
iid 2(BC)(CD) 


=> BD =./BC? + CD2 — 2(BC)(CD)cos6, 





Setting all the corresponding values in above expression, we get 
Figure 2: In plane quadrilateral ABCD, 
the sides AD, BC, CD & the angles 0; & 
0p are known 








5 92 : (“ “) 592, (“ “) 
= 222 2(24 "4 )_ 222) ot 2 eo 
R |Asin 5 + 4sin 5 8sin 5 sin 5 


| 
=R|2-2 (1 — 2sin2 22) ( — 2sin2 G5) 


BD = R, 2 — 2cos@,cos(Az — a4) sitet anceen CO) 





Let <“ADB=65 > <BDC= <ADC - ADB =06,-6, <BCD = 6¢. Now applying sine rule in 
plane ABCD (see figure 2 above) as follows 
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sin=<BDC  sin<BCD n(0) — 8) BC sin6¢ 
OF > = = — 
BC BD euie BD 


Setting all the values (take the value of BD from (3) ), we get 


2Rsin +4 es y1-—cos?6- 
RV2- Tee —A,) 


sin(@p —6) = 





2sin2 be 1 (sin sin (2 5 shy) 


SE  ————————————————— (seeting value of cos@, from (2) ) 


2 — 2cos@,cos(A, — A,) 


es | — gin? 22 cin2 (A2 =A 
: 2sin= _J/1—sin Pe sin2 (4254 7) ) 


2 — 2cos@,cos(A, — A,) 


6) — 6 = sin 





es 292 oin2 (42 = 41 
2sin+_/1— sin Pe sin2 (42 7) ) 


2 — 2cos@,cos(A, — A,) 


6 =6, —sin + 


Since 0 < Op < _ hence taking cosine both the sides, we get 


be | — cnr 22 cin? (42241 
- 2sin* /1—sin Po sin2 (42 7) ) 


2 — 2cos@zcos(Az — A,) 


es aes 2 2 ev) Ang —A4 
s 2sin = ja sin* “> sin (254) 


/ 2 — 2cos@,cos(A, — A,) 





cosé = cos} 0p —sin 


cosé = cos@pcos |} sin 


be 22 o92 (A2 = 
2sin = — sin Pz sin? (252 7) ) 


2 — 2cos@,cos(Az — A,) 


+ sin@psin} sin7+ 


Setting the value of cos@p from eq(1), 





2 — 2cosp2cos(A, — A,) — 4sin? mae — sin? es sin? (2 7) hy) 





2 — 2cos@,cos(A, — A,) 


in B 22.2 (42 = 
nt [sont C8) [a 
— sin a ee =e 


+ SSS 
2 — 2cos$,cos(A, — A,) 
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2—2(1—2sin? £2) (: — 2sin? (25% 4) — 4sin? 22.4 asint 8 sin? (254) 


2 — 2cosd@zcos(Az — A,) 


2sin ja — sin? ® sin? (254) 1 (* =~) 


= sin Pa sin (25% = “) 
2 2 


— sin? 5 sin? 


2 — 2cos@,cos(A, — A,) 





(254 z 41) — Asin? sin? (254 7) Ay 1) + Asin* = be sin? (254 54) 


= sin ca sin (=) 
2 2 


2 — 2cos@,cos(A, — A;) 


2sin es ja + sin* = $e — 2sin? oe o aan 
291. 4 (427% 
+ ———_—________——— |] 1 — sin 5 sin (=) 


J2- —- —A,) 


db, 2sin (2 z A) (1 — sin? bay 


CP) SS 
= = sin->-sin —— 
2 2 — 2cos@2cos(A, — A,) 


2sin $2 [1 — sin? Pe sin? (“2 54) rs ee 
291. (427% 
reat 2, aaa aes 1—sin 5 sin (=) 


2 — 2cos@,cos(A, —A,) 


2sin Bh sin? (254) |cos* & + 2sin % 














Ae +) 


(4) 





Now, using cosine rule in AADB as follows 


AD? + BD? — AB? 
2(AD)(BD) 


cos_<ADB = 


=> AB? = AD? + BD? — 2(AD)(BD)cosé 
Setting all the corresponding values in above expression, we get 


AB? = (2Rsin aye (R 2 — 2cosdzcos (A, — a) —2 (2nsin ) (R 2 — 2coszcos (A, — a) cosé 


AB* 2s #1 
Re = = Asin? — + 2 — 2cosd,cos(A, — 44) — Asin cosd 2 — 2cos¢,cos(A, —A 
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Setting the value of cosd./2 — 2cos@,cos(A, — A,) from eq(4), we get 


AB* 2 Pt 
— = Asin’? = + 2 —2cosd,cos(Az — 4) 


R2 
Ao—A 
_ asin & asin sin »(2—*) cos? & 


. 22 5 At ne (2) 9 G2 ne (25) 
+ 2sin 9 1—-—sin au 5 1-— sin 3 sin 3 


AB? | 
— = Asin? Ps eee (1 — 2sin? 2) ( — 2sin? @>) 





R2 





AB? Ao—A Ao —A Ao—A 
— = sin? fs + sin 2 fe + sin? (=| — 2sin? Pe sin? (=) — 2sin? P sin? (=) 


4R2 
2 Pt ey: ne (2) 
2 2 2 


fi. be 


— 2sin — sin — 
sin 2 sin 9) 


591 (“ =) ; 
= sine 1 — 2sin? as + sin 


+ 2sin? Pt sin? Pe sin? (+ 


pi. be 


— 2sin — sin — 
sin 2 sin 2 








A, —42 
= sin? Pt cosy —2,) + sin? B costa, —2,) + sin? ( < +) (1 + 2sin? 2 Ps sin? 22) 


pi. be 


— 2sin — sin — 
sin 2 sin 2 
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A,—A 
7 (sin? ms + sin? Pa) cos(A, —a,) + (1 + 2sin? Ps sin? &) sin? (=) 





_ Py 
-—2? — — 
sin-> sin 
AB? 
aR? = C (let any arbitrary constant) 
AB 
oR vc (5) 


Where, the value of constant C is given as 
A,—A 
C= (sin? Bs + sin? Pe) cos(A, —a,) + (1 + 2sin? Ps sin? f) sin? (=) 


Sint 2) (1 sine sine (25%) 


a 





Now, in isosceles AAOB (see figure 1 above) 


length of great circle arc AB 
OA=OB=R, <= A0B = eo? 


AB 
< AOB = 2si (=) 
sin“ |p 


Setting the value of AB/, » from eq(5), 


length of great circle arc AB 
—< AOB = song of grememdearce’’ = 2sin-1(VC) 


length of great circle arc AB = 2Rsin-1(VC) 


Since, great circle arc AB is the minimum distance between two given points A & B on the sphere hence 





Where, C is HCR’s constant or distal constant given as 





V 0<¢,,¢2,|4, -A,|<a >0<C<1 


NOTE: It is obvious that the value of constant C depends on the difference of angles of longitude A, — A, 
rather than the individual values of A, & A, measured from a reference plane (like prime meridian for the 
globe) hence if the difference of angles of longitude is AA then setting 2, — A, = Adin the expression of distal 


constant C, we get 
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V O<¢,,¢,AA<n >0<CK<1 


NOTE: It’s worth noticing that the above formula of HCR’s constant C has symmetrical terms i.e. if 6, & dz are 
interchanged, the formula remains unchanged & hence the value of C is unchanged. It also implies that if the 
locations of two points for given values of latitude & longitude are interchanged, the distance between them 
does not change at all. 


Since the equator plane divides the sphere into two equal hemispheres hence the above formula is applicable 
to find out the minimum distance between any two arbitrary points lying on any of two hemispheres. So for 
the convenience, the equator plane of the sphere should be taken in such a way that the given points lie on 
one of the two hemispheres resulting from division of sphere by the reference equator plane. 


T 
Since the maximum value of sin1(VC ) is 3 hence the max. of min. distance between two points on a sphere is 


Tt 
=2R (5) = 1R = half of the perimeter of a great circle passing through given points 


Case 1: If both the given points lie on the equator of the sphere then substituting @; = 2 = 0, we get 


Aa 
C = (sin?(0) + sin?(0))cosAd + (1 + 2sin?(0)sin?(0))sin? So 





AA AA 
— 2sin(0)sin(0) (1 — sin? (0) sin? =) (1 — sin? (0)sin? =) 


_AAY AA (si ox") 
maaan a ee since, 5 = > 





Aa 
C=sin?— >vVC= 
2 
Hence, the minimum distance between the points lying on the equator of the sphere of radius R 
AA Aa 
= 2Rsin-1(VC) = 2Rsin“1 (sin =) = 2R>-= RMA (0<AA<z) 


The above result shows that the minimum distance between the points lying on the equator of the sphere 
depends only on the difference of longitudes of two given points & the radius of the sphere. This can easily be 
proved by using simple geometry. 


If both the given points lie diametrically opposite on the equator of the sphere then substituting AA = zz in 
above expression, the minimum distance between such points 


RAA = R(t) = 7R = half of the perimeter of a great circle passing through thegiven points 


Case 2: If both the given points lie on a great circle arc normal to the equator of the sphere then substituting 
AA = 0 in the formula of distal constant C, we get 
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cS (sin? Bs + sin? Pa) cos(0) + (1 + 2sin? Ps sin? & sin? (0) 





- sin at Pe (1 — sin? Pe sin? 0) (1 — sin? Pa cin? (0) 
2 2 2 2 
= sin? os + sin? Pa _ 2sin 91 cin $e = (sin —sin -. 
2 2 2 2 2 2 
=VC= sin © — sin © 








Hence, the minimum distance between two points lying on a great circle arc normal to the equator of the 
sphere of radius R 


91. dp 
UID ge ag 


= 2Rsin-1(VC) = 2Rsin“} 5 


(0 < $,,¢2 <7) 








Illustrative Example 


Consider any two arbitrary points A & B having respective angles of latitude @, = 40° & @, = 75° & the 
difference of angles of longitude AA = 55° ona sphere of radius 25 cm. Now substituting the corresponding 
values, the distal constant C is given as follows 




















40° 5° 40° 75° 55° 

¢ = (sin? 5 + sin? 5 ) coss5° + (1+ 2sin” 5 sin? 5 ) sin? 5 
ste ee 
sin— sin 





= 0.279657662 + 0.231697606 — 0.394612375 = 0.116742892 


Hence, the minimum distance between the given points A & B 
= 2Rsin-1(VC) = 2 x 25sin-1(V0.116742892) ~ 17.43501363 cm 


The above result also shows that the points A & B divide the perimeter = 27(25) ~ 157.0796327cm of the 
great circle in two great circles arcs (one is minor arc AB of length ~ 17.43501363 cm & other is major arc AB 


of length = 139.6446191 cm) into a ratio ~ as Freee ret =1:8 


Conclusion: It can be concluded that this formula gives the correct values of the great circle distance because 
there is no approximation in the formula. This is an analytic formula to compute the minimum distance 
between any two arbitrary points on a sphere which is equally applicable in global positioning system. This 
formula is extremely useful to calculate the geographical distance between any two points on the globe for the 
given latitudes & longitudes. This is a highly precision formula which gives the correct values for all the 
distances on the tiny sphere as well as the large sphere like giant planet if the calculations are made precisely. 


Note: Above articles had been derived & illustrated by Mr H.C. Rajpoot (B Tech, Mechanical Engineering) 
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